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THE PATH OF LIGHT IN A GRAVITATIONAL FIELD. 

By H. S. UHLER, Yale University. 

While reading the recent books on the generalized principle of relativity, 
Einstein's theory, etc., I failed to find any numerical data which would give, even 
approximately, a quantitative idea of the angular deviation at various points 
on the curved ray which represents the course taken by a light wave in coming 
from a remote star to the earth after having passed close to the limb of the sun. 
All of the books examined give the total deviation, 1 im/R = 1.745", from infinity 
to infinity (that is, the angle between the asymptotes to the ray) but nothing 
more about numerical angles. Since my investigation of this problem may 
contain something of interest to the readers of the Monthly, it seems appropriate 
to present the details of the analysis in this place. 

As a purely mathematical problem, we shall first find the polar equation of 
a ray of light in an unlimited medium in which the index of refraction is a linear 
function of the reciprocal of the distance of the wave from the pole or center of 
symmetry. We shall then simplify the equation by adapting it to the special 
case of a gravitational field. Finally, we shall obtain so much of a series expansion 
as will be necessary for the calculation of the very small angular deviations in- 
volved. 

In order to avoid undue distraction in the midst of the subsequent analysis 
we shall now state a theorem, 2 borrowed from geometrical optics, which will 
constitute the keynote to the plan of attack to be followed in this paper. For a 
transparent medium having the property that the index of refraction possesses 
spherical symmetry with respect to a given center, the product of the index of 
refraction at any point on a ray of light in the medium by the length of the 
perpendicular dropped from the center of symmetry upon the tangent to the 
ray at the point in question is constant. Let n\, th, n 3) • • • denote the indices 
of refraction at points along the same ray and let p\, p 2 , Pz, • • • symbolize the 
lengths of the corresponding perpendiculars, then 

nipi = nzpi = nzpz = • • • = constant. (1) 

This fact is independent of the functional relation between the index of refraction 
and the distance from the center of symmetry. 

In the figure let S denote the center of symmetry and let ABC represent the 
curved ray of light. B is the point on the ray nearest to S. Take the line SB 
as the axis of polar coordinates (r, 0) of any point P on the ray. SP = r, 
Z BSP = 0. Draw tangents BT and FP to the curved ray at the points B 

1 B denotes the radius of the sun and m a quantity which Eddington calls " the mass of the 
sun in astronomical units," while Weyl styles it the " gravitational radius." m is always expressed 
in terms of a unit of length, usually the kilometer. 

2 Cf. R. S. Heath, Geometrical Optics, 1887, p. 328. 
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and P, respectively. Let the lengths of the perpendiculars (SB and SF) dropped 

from S upon these tangent lines be denoted 

^ T by po and p. Let <j> represent the acute angle 

which the radius vector to the point P makes 

'•• <£ ^^T r w '*^ the normal to the ray at P. Then 
■ ^^/P^vc z Fsp = ^ The dev j ation of the ray> D> 

at the point P will be defined as the angle 
which the tangent FP makes with the apsidal 
tangent BT. D = / TIP = £ BSF = 6 - <f>. 
In general 




o v 



1 dr 
tan <j> = 

From the figure 



tan* = -.^- (2) 



FP Vr 2 - p 2 
tan d> =-=-= *- . (3) 

Elimination of <£ from relations (2) and (3) gives 

In order to integrate equation (4), p must be given as a function of r. As 
stated in the second paragraph let 

n=h + y (5) 

For the point B we may write n = A + £/p 0> Hence, according to lemma 
(1) we have p = r(hp + l)/(hr + I). 

Substituting the last expression for p in equation (4) we obtain 

d0= dr _ • (6) 

Putting u = 1/r and Wo = 1/po in equation (6) we find 

(Jh)uod6 _ - du 



VA + 2lu L _ \ (h + 2lu p)u - WJ 

V L u (h+lu ) J 

Substituting z for the expression between the brackets, the last equation 
reduces to the simple form 

■ylhjh + 21uq) _ -da , 

h + lu a Vl - 2 2 ' 

whence 

Mh + 2§T) 
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When 6 = 0; r = p , u = «o, and z = 1, therefore c = 0. Consequently the 
equation of the curved ray may be written as 

1 = n \ 9A (' + (**»+ cos \^S+Me]\ • (7) 

We shall next apply equation (7) to the special case of the gravitational field 
associated with the sun. The point S of the figure may now be considered as 
situated at the center of this star. 

According to Einstein's theory of gravitation the index of refraction at the 
point P is given by n = 1 + 2m/r to the first order of approximation in m/r. 
As a matter of fact, m^/r 2 would be negligible in the case of the sun since, for this 
mass, m = 1.473 km.,p = 6964 X 10 2 km. (corresponding to the sun's limb) and 
m 2 /p 2 = 4.474 X 10~ n . Consequently we may obtain an approximate equation 
of the path of a wave of light in the sun's gravitational field by neglecting all 
powers of m/r higher than the first in equation (7), and by putting h = 1 and 
I = 2m. [Vide relation (5).] 

Under these conditions equation (7) reduces to 

1 2m , go - 2m ,. 

- = — 5 H 5 — cos 0. (8) 

r p<f Po 2 

This is the familiar equation of one branch of an hyperbola referred to the 
nearer focus as pole. The eccentricity of the hyperbola is given by e = po/2m — 1 
which, in the practical case under consideration, has the very large value 
2364 X 10 2 .- The transverse semi-axis = (2mp )/(po — 4m) = 2.946 km., approxi- 
mately. 

For an infinite value of r equation (8) shows that 6 is obtuse and that the 
approximate value of the total angular deviation, {2B X ), from asymptote to 
asymptote, equals 4m/po- (Vide supra.) 

In order to calculate the deviations at various points along the curved ray 
it will be necessary to derive a special formula, since the numerical values of the 
deviations are so extremely small as to preclude the use of equation (8) in con- 
junction with ordinary logarithmic-trigonometric tables. 

This may be accomplished in the following manner. Let 

— P — Po + ^m __ po 2 — 2mr _ 

— r ~ r + 2m — r(p — 2m) 

By referring to the figure and equation (8) we readily see that 
D = — = cos -1 b — cos -1 a, 

or D — sin -1 a — sin -1 b, 

or D = sin" 1 [a Vl - V- - b VI - a 2 ]- (9) 

To the first order in m/po and m/r we obtain 
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where 

po + 2m 2mpo 2 

c = 



r rV + Po) 

Hence, relation (9) reduces to 



n 2m / /p \ 2 1296 X 10 3 m ,- " 

D = — \/l - I — J = V(r - p )(r + po) • 



(10) 



The first form of relation (10) shows that the deviation very rapidly approaches 
the asymptotic limit 2m/po as r increases, and the second form is especially 
convenient for logarithmic computation. 

In calculating the data collected in the table given below the values 1.473 
km. and 6964 X 10 2 km. were taken for m and po, respectively. Thus the table 
pertains to a ray that has grazed the sun's limb. For integral values, was 
taken as the independent variable and r was calculated from equation (8). The 
numbers given for r in the table have been rounded off. The angles tabulated in 
the right hand column were obtained by the aid of formula (10). 

6° r km. D" 

6964 X 10 2 

5 6991 X 10 2 0.07605 

10 7071 X 10 2 0.15152 

15 7210 X 10 2 0.22584 

20 7411 X 10 2 0.29841 

25 7684 X 10 2 0.36876 

30 8041 X 10 2 0.43628 

35 8501 X 10 2 0.50048 

40 9091 X 10 2 0.56087 

45 9849 X 10 2 0.61700 

50 1083 X 10 3 0.66842 

55 1214 X 10 3 0.71476 

60 1393 X 10 3 0.75566 

65 1648 X 10 3 0.79081 

70 2036 X 10 3 0.81994 

75 2691 X 10 3 0.84283 

80 4010 X 10 3 0.85931 

85 7990 X 10 3 0.86925 

On 5794 X 10 4 0.87250 

B E 1495 X 10 5 0.87256 

90 1646 X 10 8 0.87257 

0« °o 0.87257 

The fourth and third rows from the bottom of the table pertain respectively 
to the orbits of the planets Mercury and Earth. d M = 89° 18' 41.5", 6 E = 89° 
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44' 0", and $«>= 90° + D». The data in the last four rows fulfill the expectation 
that the curved ray has sensibly reached its asymptote not only for the earth's 
orbit but even for an observer at the distance of the inferior planet. 1 



GRAPHICAL DISCUSSION OF THE ROOTS OF A QUARTIC 

EQUATION. 

By E. L. REES, University of Kentucky. 

It is the purpose of this note to give a graphical study of the conditions which 
determine the nature of the roots of a quartic equation. Using the reduced 
form f(x) = x 4 + qa? + rx + s = 0, with q, r and s real and with discriminant 
A, we have types of the quartic for which the following are criteria 2 regarding the 
nature of the roots: 

A < 0, roots distinct, two real, two imaginary; 

A > 0, roots distinct, all real or all imaginary; 

q 2 
q < 0, s > j , roots imaginary; 

s < 7 , roots real; 

q = 0, roots imaginary; 
A = 0, at least two equal roots; 

q 2 
q < 0, s > -. , two equal real roots, two imaginary; 

— -zrz < s < j , roots real, two and only two equal; 

* — r > two pairs of equal real roots; 

s — — r^ , roots real, three equal; 

q > 0, s > 0, r 4= 0, two equal real roots, two imaginary; 

q 2 . . 

s — — r = 0, two pairs of equal imaginary roots; 

s = 0, two equal real roots, two imaginary; 
q = 0, s > 0, two equal real roots, two imaginary; 

* = 0, four equal real roots. 

The discriminant is the product of the squares of the differences of the roots 

1 The deviation of a ray from its original direction at the star is obtained by adding 2m/po to 
the value of D given in the table, and the total deviation, 1.745", from asymptote to asymptote 
is twice the value of D given in the last two lines of the table. — Editors. 

! Compare L. E. Dickson, Elementary Theory of Equations, 1914, p. 45. 



